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1 . INTRODUCTION 


The  purpose  of  this  research  study  has  been  to  devise  a physi- 
cally  sound  alternative  to  the  Smith  - Van  Ingen4  e9  procedure 
for  predicting  boundary-layer  transition.  The  need  for  develop- 
ing an  alternative  procedure  stems  from  the  following:  On  the 
one  hand,  there  is  little  doubt  that  stability  theory's  Tollmien- 
Schlichting  waves  exist  and  play  an  important  role  in  the  initial 
stages  of  transition.  On  the  other  hand,  because  the  end  state 
of  the  transition  process  is  a (highly  nonlinear)  turbulent  flow, 
linear-stability  theory  must  break  down  at  some  point  between 
that  of  the  initiation  of  Tollmien-Schliehting  waves  and  the 
transition  point  (defined,  for  example,  as  the  point  where  skin 
friction  achieves  a minimum).  In  other  words,  linear-stability 
theory  is  inapplicable  in  the  post-critical  stages  of  transition 
and  therefore  has  no  natural  way  of  specifying  the  actual  transi- 
tion point.  The  empirical  e9  method  ignores  this  conceptual 
limitation  and  thus  has  a questionable  physical  foundation. 


Under  sponsorship  of  the  Office  of  Naval  Research  (ONR) , DCW 
Industries  has  taken  an  important  first  step  toward  developing  a 
physically-sound  alternative  to  the  e9  method  for  predicting 
transition.  The  basic  approach  has  been  to  devise  a synthesis 
of  linear- stability  theory  and  turbulence-model  techniques  by 
using  each  theory  only  in  the  regime  where  its  inherent  limita- 
tions determine  its  applicability.  Our  first  step  has  consisted 
of  showing  that  a key  closure  coefficient,  A,  in  the  turbulence- 
model  equations  is  relatively  insensitive  to  frequency  after  a 
boundary- layer  disturbance  has  been  amplified  to  about  e1*  times 
its  initial  value,  provided  the  boundary  layer  is  unstable. 

While  results  of  this  first  step  have  been  very  encouraging, 
analysis  has  been  confined  to  the  Blasius  boundary  layer. 


primarily  because  of  difficulties  encountered  in  acquiring  and 
debugging  modifications  to  the  Mack  stability  program.  Addi- 
tional research  has  been  needed  to  confirm  our  hypothesis  that  a 
linear-stability/ turbulence-model  synthesis  is  feasible. 
Specifically,  two  things  need  to  be  done.  First,  for  the 
constant  pressure  case,  actual  transition  predictions  must  be 
made  using  the  stability-predicted  X profile  and  appropriate 
initial  and  boundary  conditions.  Second,  computation  of  this 
key  closure  coefficient  when  the  pressure  varies  is  needed  to 
provide  a more  stringent  test  of  our  hypothesis.  Both  ends  have 
been  accomplished  in  this  research  study. 

Section  2 reviews  the  basic  formulation  underlying  our  theoret- 
ical approach.  In  Section  3,  the  Blasius  boundary  layer  is 
analyzed  including  comparison  of  predicted  and  experimentally- 
measured  transition  points.  Effects  of  both  adverse  and  favorabl 
pressure  gradient  are  addressed  in  Section  4.  Results  and  con- 
clusions follow  in  Section  5. 


2.  THEORETICAL  FORMULATION 


In  this  section,  for  the  sake  of  completeness,  we  first  review 

4 

the  overall  approach  devised  by  Wilcox.  Then  the  turbulence 
model  equations  underlying  the  approach  are  presented.  Finally, 
the  manner  in  which  turbulence-model  parameters  are  computed  from 
a given  linear-stability  solution  is  delineated. 

2.1  OVERVIEW 

Conceptually,  our  overall  approach  to  predicting  boundary-layer 
transition  consists  of  three  interrelated  phases.  In  the  first 
phase  conventional  linear-stability  computations  are  performed  up 
to  and  beyond  the  point  at  which  a disturbance  becomes  unstable. 

As  a key  feature  of  these  first-phase  computations  we  have  found 

4 

empirically  that  the  computations  are  most  appropriately  termi- 
nated when  the  initial  disturbance  has  been  amplified  by  a factor 
of  e4.  In  the  second  phase  the  linear-stability  solutions  at  the 
e4 -amplification  point  are  used  to  compute  key  turbulence  model 
parameters  which  are  needed  to  specify  the  closure  coefficient  X 
as  well  as  initial  and  boundary  conditions  for  a subsequent 
turbulence-model  calculation.  In  the  third  phase  a boundary- layer 
program  incorporating  the  Wilcox-Traci  turbulence-model 
equations  is  used  to  predict  the  actual  transition  point.  Most 
notably,  a range  of  frequencies  is  considered  in  all  phases  of 
the  computational  procedure. 

The  rationale  of  this  approach  is  two-fold.  First , using  linear- 
stability  theory  in  the  initial  phase  of  the  procedure  allows  us 
to  take  proper  account  of  the  transition  phenomenon's  sensitivity 
to  the  frequency  of  the  disturbance.  Using  the  turbulence-model 
equations  is  quite  inappropriate  in  the  initial  phase  as  long- 
time averaging  (which  is  an  integral  aspect  of  the  turbulence- 
model  theory)  precludes  explicit  representation  of  frequency 
effects.  Second,  using  the  turbulence-model  equations  in  the 
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final  phase  allows  us  (at  least  on  a conceptual  level)  to 
represent  nonlinear  processes  which  ultimately  must  become 
dominant  as  the  boundary  layer  nears  transition  to  a turbulent 
state.  Using  linear-stability  theory  is  quite  inappropriate  in 
the  final  phase  as  the  linearization  precludes  explicit  repre- 
L sentation  of  nonlinear  effects. 

The  astute  reader  will  notice  that  a key  premise  is  implicit  in 
the  rationale  of  our  approach.  Specifically,  for  this  approach 
to  work,  it  must  be  the  case  that  either  (a)  the  transition 
phenomenon  is  frequency  independent  beyond  the  e ‘‘-amplification 
point  or  (b)  a method  exists  whereby  frequency  effects  can  be 
explicitly  expressed  in  the  turbulence-model  equations  during 
the  latter  phases  of  transition.  As  will  be  demonstrated  in 
this  report,  while  the  transition  process  is  somewhat  frequency 
dependent  beyond  the  e4-amplif ication  point,  this  frequency 
dependence  can  indeed  be  explicitly  expressed  in  the  model 
equations . 

The  remainder  of  this  section  is  devoted  to  a description  of  the 
turbulence-model  equations  and  the  manner  in  which  key  turbulence- 
model  parameters  are  computed  from  a linear-stability  solution. 

2.2  TURBULENCE-MODEL  EQUATIONS 

The  two-equation  model  of  turbulence  developed  by  Wilcox  and 
Traci^’6  is  used  in  this  study  for  predicting  boundary- layer 
development  during  the  latter  phases  of  transition.  For  incom- 
pressible boundary  layers  this  model  consists  of  the  long-time 
averaged  conservation  of  mass  and  momentum  equations  and  two 
additional  rate  equations.  Denoting  arc  length  and  surface- 
normal  distance  by  x and  y with  corresponding  velocity  components 
u and  v,  the  four  equations  of  motion  for  incompressible  boundary 
layers  are: 
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where  p is  pressure  and  v is  kinematic  viscosity.  The  quantity  e 
is  kinematic  eddy  viscosity  which  is  computed  as  the  ratio  of 
specific  turbulent  mixing  energy,  e,  to  turbulent  dissipation 
rate,  u,  i . e . , 

e = e/u  (5) 

Also,  l is  the  turbulence  length  scale  computed  according  to 

l = eVw  (6) 

Finally,  the  quantities  a,  a*,  B,  B*,  o and  a*  are  closure 
coefficients  whose  values  appropriate  for  turbulent  flows  have 
been  found  empirically^  to  be  given  by  the  following: 

I B = 3/20,  B*  = 9/100 

a = 1/2,  a*  = 1/2 

a - j [l  - (1- X)  exp  (-ReT/2)] 

] a*  [i  - Cl-X)  exp  C-2ReT)] 

Where  Re^  is  Reynolds  number  based  on  turbulence  properties 
defined  as 

r 

• 5 
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ReT  = e'H/v  (8) 

and  the  quantity  X is  an  additional  closure  coefficient  whose 
value  has  been  found  by  Wilcox^  to  vary  with  freestream 
turbulence  level,  surface  temperature  and  pressure  gradient. 

For  the  Blasius  boundary  layer  a value  of  X = 1/11  has  been  found 
to  yield  transition-point  predictions  in  reasonably  close  agree- 
ment with  experimentally  measured  transition  points. 

Experience  with  the  model  equations^  has  shown  that  the  model's 
ability  to  predict  transition  hinges  most  heavily  upon  the 
closure  coefficient  X. 

Equations  (l)-(4)  must  be  solved  subject  to  boundary  conditions 
at  the  surface,  y=0 , and  at  the  boundary- layer  edge,  y=5,  where  6 
is  boundary- layer  thickness.  Prior  studies  have  established  that 
for  perfectly-smooth  surfaces  (rough  surfaces  are  not  considered 
in  this  study) , we  must  require 


u = v = 

e = 0 

at 

O 

II 

X 

(9) 

and 

^ 20 
“ - T 

V 

y2 

as 

y -*■  0 

(10) 

Denoting  conditions  at  the 

boundary- layer 

edge  by 

subscript 

we  must  also 

require 

u = Ue  (x) 

at 

y = 6 

(ID 

while  e fx) 
e 

and  a)  (x)  must 
e 

satisfy  the 

following 

simplified 

equations : 

de 

U ® 

e 35 ~ 

= - e 

e e 

> at 

y = 5 

(12) 

do 

U e 

- - go)  3 1 

e 35 T 

e 1 

2.3  COMPUTING  TURBULENCE -MODEL  PARAMETERS 
FROM  A LINEAR-STABILITY  SOLUTION 


As  stated  in  Subsection  2.1,  before  the  turbulence-model  phase 
of  the  computation  can  be  initiated,  initial  and  boundary  con- 
ditions as  well  as  the  closure  coefficient  A must  be  deduced  from 
the  linear-stability  solution.  To  do  so  we  must  first  note 
several  important  facts  about  the  various  turbulence-model 
parameters . 


First  of  all,  in  order  to  define  initial  and  boundary  values  for 
e and  oj  we  note  that  Wilcox  and  Chambers  have  defined  these 
quantities  in  terms  of  the  fluctuating  vertical  velocity,  v', 
according  to  the  following: 

e = | <v  * 2 > (13) 

and 

T.l  s ( o.r  * / 2 ^ 

(14) 


0) 


. 3v  <(3v'/3y)  > 
g*  <v'z> 


where  < > denotes  long-time  average.  As  the  linear-stability 
solution  yields  the  instantaneous  v'  profile,  Equations  (13)  and 
(14)  are  sufficient  to  determine  e and  u>  profiles  provided  a 
suitable  time-averaging  process  can  be  defined. 


Next,  in  order  to  determine  A,  we  first  rewrite  Equation  (3) 
under  the  standard  parallel-flow  assumption  (i.e.,  vso)  and  in 
the  limit  of  small  turbulent  Reynolds  number  (i.e.,  ReT  <<  1) . 
The  result  is 


u 


3e  _ 3 . 

3x  ITT 


3u 

3y 


e - 


S*we  + 


(15) 


+ Note  that  in  the  actual  turbulence-model  phase  of  the  compu- 
tation, Equation  (3)  is  used.  Equation  (15)  has  been  introduced 
only  to  devise  a method  for  computing  A from  the  linear-stability 
solution. 
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Now,  noting  the  definition  of  e given  in  Equation  (13)  , we  can 
derive  an  exact  equation  for  its  evolution  by  taking  the  v' 
moment  of  the  v-momentum  equation  and  time  averaging.  Assuming 
that  in  the  limit  Re^.-*-0  the  triple  correlation  term  is  negligibly 
small,  we  obtain 


<v'l£ 

p 3y 


> 


(16) 


Comparison  of  Equations  (15)  and  (16)  shows  that  the  closure 
coefficient  A can  be  defined  as  follows: 


\ 
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<v  ' 2 > 


2 


v»  3 p ' 
p Zy 


I 3u/3y  I 


> 


(17) 


The  final  step  in  relating  the  linear- stability  and  the  turbulence- 
model  parameters  is  to  define  an  appropriate  long-time  average. 

To  do  so  we  note  first  that  in  the  linear- stability  solution  the 
velocity  and  pressure  fluctuations,  v'  and  p',  are  written  as 


v *(x,y,z,t) 
p *(x,y,z,t) 


= 4>  (y ) exp 

= (y)  exp 


AAA 

ax+gz-oo 


t) 


A /V  A.  ) 

ax  + 3z  -cot , 


(18) 

(19) 


where  t denotes  time,  z is  distance  normal  to  the  x and  y axes, 

A A 

is  freestream  velocity,  p is  density,  a and  3 are  wave  numbers, 
to  is  frequency,  and  the  functions  $ (y)  and  tt  (y)  are  the  complex 
amplitude  functions  of  the  fluctuating  flow  variables  v'  and  p'. 


To  evaluate  the  time-averaged  quantities  appearing  in  Equations 
(13),  (14)  and  (17)  we  use  the  following  deiinition: 

t+T 

(x,y,z,T)dt  (20) 

-T 


<y>  = lim 


(oT-*-«>  2T 


t + 

M*' 

Vt-i 
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Then,  working  with  the  real  parts  of  the  linear-stability 
solution,  performing  all  time  averages  indicated  in  Equations 
(13),  (14)  and  (17),  and  denoting  local  Reynolds  number  by  R 
yields  the  following: 


e/U^ 

■ l{*l*  *l) 

(21) 

(d<j>r/dn) 

voj  3 

2 + (d<j>i/dn)2 

(22) 

u2  b#r2 

i * *i) 

R_1  [(d4>r/dn)2  + 

(d^/dn)2]  - 2 1 

<j)rdTTr/dn+4>idTTi/dnJ 

(*r  ♦ *i)  1 

3U/9n|  (23) 

where  U = u/U^,  subscripts  r and  i denote  real  and  imaginary  parts, 
and  q is  dimensionless  distance  defined  by 


(24) 


with  Ug  denoting  horizontal  velocity  at  the  boundary- layer  edge. 
Note  that  for  the  constant-pressure  case,  Ug  * U^.  Equations 
(21) -(23)  are  the  desired  relations  which  can  be  used  to  define 
X and  the  initial  and  boundary  conditions  for  a turbulence-model 
computation  in  terms  of  a given  linear-stability  solution. 
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3.  RESULTS  FOR  CONSTANT  PRESSURE 


In  this  section  we  implement  the  basic  method  for  the  constant- 
pressure  case,  i.e.,  for  the  Blasius  boundary  layer.  First, 
we  compute  the  closure  coefficient  X for  a range  of  frequencies 
and  demonstrate  that  beyond  the  e4 -amplification  point  X 
asymptotes  to  a nearly  frequency-independent  profile.  Next, 
we  compute  e and  oj  profiles  at  the  e4 -amplification  point 
which  can  be  used  to  initiate  a turbulence-model  computation. 
Then,  we  establish  boundary  conditions  for  the  turbulence- 
model  phase.  Finally,  transition  predictions  are  made  and 
compared  with  experimental  data, 

3.1  EVALUATION  OF  THE  CLOSURE  COEFFICIENT  X 

4 

Using  Equation  (17),  Wilcox  has  computed  X profiles  for  the 
Blasius  boundary  layer.  For  the  sake  of  completeness  results 
obtained  in  the  Wilcox  study  are  repeated  in  this  subsection. 

A large  number  of  linear-stability  computations  have  been  per- 
formed  with  the  Mack  stability  program;  all  computations  have 
been  done  with  the  spatial  amplification  theory  option.  Both 
Reynolds  number  and  frequency  have  been  varied  in  order  to 
determine  the  variation  of  X throughout  the  Reynolds -number/ 
frequency  plane. 

Figure  1 shows  computed  X profiles  at  nine  Reynolds  numbers 
corresponding  to  one  stable  case,  one  neutrally-stable  case, 
and  seven  unstable  cases  corresponding  to  amplification  from 
the  neutral  case  by  factors  of  e11  with  values  of  -n  ranging  0 
to  10;  for  all  nine  cases  the  frequency  is  given  by 

Fr  = wv/U*  = 2 • 10" 5 (25) 

We  are  thus  following  the  evolution  of  X for  a constant- 
frequency  disturbance  at  a Reynolds  number  upstream  of  the 
neutral  point  corresponding  to  the  frequency  given  in 
Equation  (25)  . 
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As  shown,  although  A varies  rapidly  with  y near  the  surface,  all 
nine  curves  display  approximately  the  same  variation.  However, 
above  a value  of  n = 2,  the  various  A profiles  vary  rapidly  with  n 
and  do  so  in  dissimilar  manner  for  amplification  rates  up  to  e4'56. 
For  example,  at  the  lowest  Reynolds  number  for  which  the  solution 
is  stable,  A is  negative  above  n=  2 . As  we  move  to  the  neutral 
point  we  find  that  A vanishes  for  values  of  n in  excess  of  2.5. 

Then  as  Reynolds  number  increases,  A varies  more  and  more  rapidly 
with  n and  asymptotes  to  a single  curve  for  amplification  ratios 
in  excess  of  between  e3  and  e4. 

Figure  2 shows  similar  curves  for  a dimensionless  frequency  Fr 
given  by 

Fr  = 3-10'5  (26) 

Again,  the  curves  collapse  to  a single  curve  for  amplification 
ratios  in  excess  of  e4. 

Computations  have  been  performed  for  frequencies  covering  the 
entire  stability  diagram.  For  each  frequency  considered  the 
computed  A profiles  always  asymptote  to  a universal  profile 
beyond  the  e4-amplification  point  with  a subtle  qualification. 

That  is,  referring  to  Figure  3,  as  Reynolds  number  increases  the 
upper  branch  of  the  stability  diagram  eventually  is  reached  and 
we  again  enter  a stable  region.  As  we  approach  this  upper 
neutral  point,  the  A proflies  begin  to  fall  back  to  those  typical 
of  low  Reynolds  numbers. 

The  rapid  variation  of  A near  n=0  results  from  a breakdown  in 
the  basic  closure  approximations  near  the  surface.  That  is,  the 
the  production  term  in  the  <v' 2>  equation,  goes  to 

zero  quadratically  with  distance  from  the  surface  so  that,  in 
terms  of  n» 

<-  §£’>  ^ n2  as  n-*-0  (27) 
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Figure  2.  Computed  X profiles  for  dimensionless  frequency 
Fr  = 3*10“5;  zero  pressure  gradient. 


By  contrast,  the  modeled  production  term  for  ReT-*-Q  behaves  as 


^ x|H  e * Xn“  . as  n-0  C28) 

Consequently,  close  to  the  surface  we  have 

X -v  n 2 as  n-*-0  (29) 

This  modeling  shortcoming  is  of  little  consequence  as  dissipation 
exceeds  production  near  q = 0.  Hence,  for  the  remainder  of  this 
discussion  our  focus  will  be  upon  the  region  between  n=l  and  the 
outer  edge  of  the  boundary  layer,  n=5. 

Figure  4 shows  computed  X profiles  for  several  frequencies  and 
amplification  ratios.  As  shown,  all  of  the  computed  X profiles 
cluster  about  the  approximate  profile  defined  by 

X = .0093  + .00156  exp  j^-  (n*l)j  (30) 

The  fact  that  the  computed  X profiles  correlate  with  the  profile 
defined  in  Equation  (30)  independent  of  frequency  is  extremely 
important.  This  is  our  first  indication  that  using  the  turbulence- 
model  equations  to  describe  the  latter  stages  of  transition  is 
feasible . 

It  is  interesting  to  examine  the  rate  of  approach  to  the  asymptotic 
profile  which  is  most  conveniently  done  in  terms  of  the  average 
value  of  X defined  as  follows: 

5 

X dn  (31) 

Figure  5 shows  X as  a function  of  displacement  thickness  Reynolds 
number.  Re.*;  note  that  Re.,  denotes  the  neutral-stability  value 

0 0 w j 

of  Re^*  for  a given  frequency.  As  shown,  for  the  higher 


Figure  4.  Profiles  of  the  closure  coefficient  X for  various 

frequencies  and  amplification  ratios;  zero  pressure 
gradient . 


Figure  5 . Variation  of  X with  Reynolds  number  for  several 
frequencies;  zero  pressure  gradient. 
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frequencies  X approaches  its  asymptotic  value  most  rapidly. 

For  Fr  = 0.5‘10~5,  the  lowest  frequency  at  which  computations 
have  been  done,  the  approach  to  the  asymptotic  values  lie 
between  X=.066  and  X=.083,  as  compared  to  the  postulated^ 
turbulence-model  value  of  .091. 

Further  examination  of  the  X variation  with  Reynolds  number 
shows  that,  for  the  higher  frequencies,  the  peak  value  is 
achieved  at  an  amplification  ratio  of  about  e4  while,  for  the 
lowest  frequency,  X is  about  half  its  asymptotic  value  at  this 
ratio.  Figure  6 presents  the  variation  of  X with  frequency 
for  amplification  ratio  e4,  including  a correlation  of  the 
computed  values;  the  correlation  is: 

Tn=4  ’ TS  i1  - exP  [-  “ (104Fr)]  | (32) 

The  asymptotic  values  of  X,  denoted  as  X , are  also  shown 
for  reference.  Two  key  conclusions  can  be  drawn  from  the 
computed  variation  of  X.  First , beyond  an  amplification  ratio 
of  e4,  X varies  slowly  with  Reynolds  number.  Second,  X is  only 
weakly  frequency  dependent  for  frequencies  in  excess  of 
Fr  = 10  5.  These  two  points  lend  further  credence  to  our 
hypothesis  that  turbulence-model  equations  potentially  can 
yield  a sound  physical  description  of  the  latter  stages  of 
transition. 

3.2  INITIAL  PROFILES 

Having  established  Equation  (30)  as  a satisfactory  correlation 
of  computed  X profiles  we  now  turn  to  initial  profiles  for  e 
and  w which  are  needed  to  initiate  a boundary-layer  computa- 
tion beyond  the  e4 -amplification  point.  Figures  7 and  8 show 
computed  profiles  for  e and  ui  at  the  e4-amplif ication  point. 
These  profiles  have  been  obtained  from  the  linear-stability 
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Figure  6.  Correlation  of  the  average  value  of  X with  frequency 
for  an  amplification  ratio  of  ek ; zero  pressure 
gradient . 


Figure  7.  Computed  turbulent  mixing-energy  profiles  after  an 
initial  disturbance  has  been  amplified  by  a factor 
of  e1*:  zero  pressure  gradient. 
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computations  and  Equations  (21)-(_22).  In  each  of  the  figures 
profiles  are  shown  for  five  values  of  the  dimensionless 
frequency,  Fr. 

The  e profiles  are  displayed  in  terms  of  e/eg  vs.  y/6  where 
eg  is  the  value  of  e at  the  boundary-layer  edge  and  6 is 
boundary- layer  thickness.  The  u)  proflies  are  parameterized 
as  a function  of  4>  vs.  y/<5  where  $ is  defined  by 


* 


_ 0 — 2 
TV  n 


X(D 


(33) 


where  the  value  of  6 is  given  in  Equation  (7) . As  shown,  the 
profiles  are  similar  for  the  various  frequencies  although 
significant  differences  in  amplitude  are  clearly  indicated. 
These  e and  a>  profiles  have  been  fitted  with  cubic  splines 
which  yield  a high  degree  of  accuracy  for  values  and  slopes 
of  the  profiles. 


3.3  BOUNDARY  CONDITIONS 

Before  proceeding  to  the  turbulence-model  computations, 
boundary-conditions  must  be  specified  at  the  boundary- layer 
edge,  y=<5.  As  in  all  previous  computations,  the  boundary 
condition  for  e is  determined  by  the  freestream  turbulence 
intensity  according  to 


where  T'  is  turbulence  intensity  given  in  percent.  To  set  the 
value  of  o)e,  we  have  analyzed  the  linear-stability-predicted 
values.  Figure  9 shows  that  at  the  e4 -amplification  point, 
vue/Ug  correlates  with  the  dimensionless  frequency,  Fr, 
according  to  the  following  formula: 
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Figure  9.  Correlation  of  boundary-layer-edge  value  of  the 
turbulent  dissipation  rate  with  dimensionless 
frequency;  zero  pressure  gradient. 
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(35) 


In  Equation  (35),  <d  denotes  the  value  of  u>  at  the  e4- 

en  e 

amplification  point.  Furthermore,  according  to  the  turbulence 
model  equations,  we  satisfies  the  equation 


U 


e 


dw* 

e 

dx 


(36) 


so  that,  in  terms  of  the  function  4 defined  in  Equation  (33), 
u)q  should  vary  with  distance  from  the  plate  leading  edge  as 


where  $ is  the  value  of  $ at  y=6  and  xn  is  the  value  of  x at 
e '0 

the  e4 -amplif ication  point.  Figure  10  compares  the  linear- 
stability-computed  variation  of  u>  with  the  variation  given 
by  Equation  (37).  As  shown,  the  two  variations  are  reasonably 
close.  Hence,  Equations  (35)  and  (37)  have  been  used  in  the 
computations . 

Before  proceeding  to  results  of  the  computations,  it  is  in- 
structive to  pause  and  discuss  the  implication  of  the  derived 
boundary  condition  for  w.  Unlike  previous  turbulence-model 
transition-prediction  computations,  we  are  not  at  liberty  to 
arbitrarily  specify  oje . This  is  an  important  point  because 
the  boundary  condition  for  w has  been  ambiguous  in  the  past 
and  its  specification  has  provided  the  equivalent  of  an 
"adjustable  parameter"  in  the  theory.  Elimination  of  this 
adjustability  enhances  utility  of  the  method,  particularly  if 

9 

it  proves  to  be  an  accurate  predictive  tool. 
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Figure  10.  Comparison  of  linear-stability-computed  variation 
of  u>e  with  model-equation-predicted  variation; 
zero  pressure  gradient. 


3.4  TRANSITION  PREDICTIONS 


To  assess  the  accuracy  of  the  method,  a series  of  computations 
has  been  performed  for  the  Blasius  boundary  layer.  In  all 
calculations  we  have  used:  Equation  (30)  to  define  X;  cubic 
spline  fits  for  the  initial  e and  go  profiles  shown  in  Figures 
7 and  8;  and  Equations  (34),  (35)  and  (37)  to  define  boundary- 
layer  edge  conditions  for  e and  u>.  Further,  in  the  computa- 
tions, freestream  turbulence  intensity  varies  from  .011  to 

.50%  while  0.1  < lO^Fr  < 0.4.  The  upper  bound  on  Fr  has 

-4 

been  chosen  because  0.4-10  is  the  highest  frequency  for 
which  the  stability  diagram  is  sufficiently  wide  to  permit 
amplification  much  beyond  a factor  of  e4  so  that  our  procedure 
becomes  meaningless  at  higher  frequencies.  The  lower  bound 
on  Fr  has  been  chosen  as  the  critical  Reynolds  number  for 
Fr=10  3 is  well  above  the  transition  Reynolds  numbers  observed 

g 

by  Schubauer  and  Skramstad  , indicating  that  lower  frequencies 
are  relevant  only  for  much  lower  intensities  than  .01%.  Since 
we  are  comparing  our  results  with  the  Schubauer -Skramstad  data 
which  were  taken  for  intensities  in  excess  of  .02%,  this  is  a 
reasonable  lower  bound  on  Fr. 

Figure  11  compares  computed  and  measured  transition  momentum- 

thickness  Reynolds  number,  Refl  , for  four  frequencies,  viz, 

4 

10  Fr=0.1,  0.2,  0.3  and  0.4.  As  shown,  the  data  lie  mostly 

4 

near  the  curves  for  10  Fr=0.2  and  0.4.  This  result  is  quite 
encouraging  as  this  band  of  frequencies  is  comparable  to  that 
predicted  by  linear- stability  theory3  to  be  most  unstable  for 
the  range  of  transition  Reynolds  numbers  observed  by 
Schubauer  and  Skramstad. 


Figure  12  compares  computed  and  measured  transition  Reynolds 

number  based  on  distance  from  the  plate  leading  edge.  Re  ; 

4 xt 

only  the  curve  computed  with  10  Fr=0.2  is  shown  for 


the  computations,  viz,  that  following  the  predicted  variation 
of  Re  with  T'  down  to  about  T'=.03%  indicates  an  asymptote 

to  a finite  value  of  Re  as  T'+O,  similar  to  the  trend  of 

the  data.  However  below  T'=.03$,  the  predicted  Rex  then 
increases  rapidly  implying  that  Rex  may  approach  an  infinite 
value  as  T'  -*•  0.  Although  not  show?!  in  the  figure,  a similar 
trend  appears  for  the  other  frequencies. 

As  evidence  of  the  method's  ability  to  properly  treat  spectral 
effects.  Figure  13  compares  an  inferred  "most -dangerous" 
frequency  with  the  linear- stability-theory  predicted  "most- 
dangerous"  frequency  band.  The  inferred  curve  has  been  con- 
structed by  forcing  agreement  between  model -predicted 
transition  point  and  the  Schubauer-Skramstad  data.  As  shown, 
the  curve  lies  just  slightly  outside  the  linear-stability- 
theory  band. 

As  a final  note,  the  actual  transition  predictions  have  been 
found  to  be  relatively  insensitive  to  initial  profiles.  To 
make  a similar  claim  about  whether  starting  the  turbulence- 
model  computation  at  a different  point  (e.g.,  when  the 
amplification  ratio  is  e3  or  e5)  would  require  further 
research.  Some  computations  have  been  done  from  the  e5- 
amplif ication  point  and  this  appeared  to  have  little  effect 
on  the  predicted  transition  point. 
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4.  EFFECTS  OF  PRESSURE  GRADIENT 


In  order  to  investigate  effects  of  pressure  gradient,  we  now 
repeat  the  analysis  of  Section  3 using  Falkner-Skan  velocity 
profiles  for  both  adverse  and  favorable  pressure  gradients. 
First,  we  compute  X profiles  at  the  e4 -amplification  point 
and  seek  a frequency- independent  correlation.  Then  the 
boundary- layer  edge  condition  for  to  is  correlated  with 
frequency.  Results  follow. 


4.1  EVALUATION  OF  THE  CLOSURE  COEFFICIENT  X 


The  Falkner-Skan  profiles  satisfy  the  following  equation: 


€1  + f i!£  ♦ 6t 


[i  - - 

L \dn  / _ 


which  is  solved  subject  to: 


f = — = 0 at 

dn 


Si  - 1 as 
dn 


In  Equation  (38)  , f is  the  dimensionless  streamfunction  so 
that  the  velocity  u is  obtained  from 


u = U « 
e d* 


Also,  the  freestream  velocity,  Ug,  is  given  by 


U (x)  = Cx" 
e 


where  C is  a constant  and  m is  related  to  the  constant  0™ 

r o 

appearing  in  Equation  (38)  through  the  following  equation: 


(42) 


Finally,  the  scaled  coordinate  n is  defined  in  terms  of  x,y,v 
and  Ug  by 


(43) 


Note  that  for  zero  pressure  gradient  m=0  and  Equation  (43) 
differs  from  Equation  (24)  by  a factor  sfl.  For  the  sake  of 
consistency  with  our  earlier  computations,  all  results  in  this 
section  are  cast  in  terms  of  n as  defined  in  Equation  (24) 
rather  than  in  terms  of  n- 


For  symmetry,  we  consider  one  adverse  pressure  gradient  profile 
and  one  favorable  pressure  gradient  profile,  viz,  we  select 

1-0.18  , adverse  Vp 

(44) 

+0.18  , favorable  Vp 

which  correspond  to  m=-.083  and  m=+.099  for  the  adverse  and 
favorable  cases,  respectively. 

Turning  first  to  adverse  pressure  gradient,  Figure  14  shows 
computed  values  of  X for  several  frequencies  near  the  e4- 
amplif ication  point;  the  correlation  for  zero  pressure  gradient 
[Equation  (30)]  is  shown  for  reference.  As  can  be  seen,  the 
various  X profiles  show  a stronger  frequency  dependence  than 
exists  in  the  absence  of  pressure  gradient.  Furthermore,  the 
disparity  in  the  various  profiles  fails  to  diminish  as  n 
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Figure  14.  Profiles  of  the  closure  coefficient  X for  various 
frequencies  near  the  e4-amplif ication  point;  ad- 
verse pressure  gradient. 


increases.  Close  examination  of  the  computed  results  shows 

that  the  minimum  value  of  A increases  as  Fr  increases.  This 

trend  is  also  present  for  zero  pressure  gradient  (see 

Figure  4),  although  not  as  pronounced  as  for  8pS=-0.18.  This 

observation  suggests  that  a better  correlation  might  be 

obtained  by  working  with  the  ratio  of  A to  its  minimum  value, 

A - . Figure  15  shows  that  the  ratio  A/A  . does  indeed 
mm  ° min 

correlate  nicely;  the  correlation  is 

X . - -9488  ♦ .0512  exp  ll  ■ jVI ; eps  - -.18 

min  L ' 'J 

(45) 


Reexamination  of  the  zero-pressure  gradient  data  shows  that 
rewriting  Equation  (30)  in  terms  of  and  y/6  improves 
the  correlation;  the  revised  correlation  is 


A 

X . 

min 


.8564  + 


.1436  exp 


0 


(46) 


Turning  finally  to  favorable  pressure  gradient,  we  find 
similar  results  to  those  obtained  with  zero  and  adverse 
gradients,  i.e.,  is  weakly  frequency  dependent.  The 

favorable  gradient  correlation  is 


A 

A . 
min 


604  + .396  exp 


y (j  • *)’ 


=+.18 


(47) 


Figure  16  displays  the  three  inferred  profiles  defined 

in  Equations  (45) -(47). 


As  a final  comment,  examination  of  the  computed  values  of 

A . shows  that  an  excellent  correlation  of  A . can  be 
min  min 
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Figure  15.  Correlation  of  *Amln  with  adverse  pressure 
gradient . 
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devised  which  is  independent  of  pressure  gradient.  Specif- 
ically, is  found  to  vary  inversely  with  displacement - 

thickness  Reynolds  number,  Re^*,  according  to  the  following 
formula: 

Amin  = 38/Re<S*  (48) 


Figure  17  compares  computed  values  of  X . with  Equation  (48). 


4.2  BOUNDARY  CONDITIONS 

As  the  initial  profiles  appear  to  have  little  effect  on 
transition  predictions  attending  use  of  the  established  pro- 
cedure, we  now  proceed  directly  to  boundary  conditions.  Again 
we  use  Equation  (34)  to  define  the  boundary- layer  edge  value 
of  the  turbulent  mixing  energy.  Also,  Equation  (37)  remains 
valid  for  determining  we/we  although  the  validity  of 
Equation  (35)  must  be  ascertained. 


Figure  18  shows  computed  values  of  vue  /Ug2  as  a function  of 

Fr  for  favorable,  zero  and  adverse  pressure  gradients.  As  in 

7/4 

the  8CC=0  case,  we  find  that  vu  /U  2 is  proportional  to  Fr  . 

r u 6/>  6 

Furthermore,  the  following  modified  relation  correlates  all  of 
the  data  to  within  a few  percent. 


28  Fr 
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(49) 


Hence,  even  in  the  presence  of  pressure  gradient,  the  edge 
value  of  a)  is  uniquely  determined  as  a function  of  the 
frequency  of  the  disturbance. 


37 


0 ' 1000  2000  3000  4000  5000  6000 


Re<s» 


Figure  17.  Correlation  of  A . as  a function  of  displace 
ment  thickness  m n Reynolds  number,  Re**. 
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5 , DISCUSSION 


While  some  fine  tuning  of  the  procedure  is  needed,  the  results 
described  in  Sections  3 and  4 are  very  encouraging.  Our  pro- 
posed synthesis  of  the  linear-stability  and  turbulence-model 
transition-prediction  methods  appears  to  be  quite  sensible. 

On  the  one  hand,  using  linear-stability  theory  in  the  initial 
stages  of  the  transition  phenomenon  (i.e.,  up  to  the  e4- 
amplification  point)  provides  a physically-sound  procedure 
which  takes  proper  account  of  spectral  effects.  On  the  other 
hand,  we  believe  that  using  the  turbulence-model  equations  to 
describe  the  latter  stages  of  transition  with  starting 
profiles  for  e,  <d  and  X and  the  boundary- layer  edge  boundary 
condition  on  w as  predicted  by  linear-stability  theory  is 
more  physically  sound  than  continuing  beyond  the  e4 -amplification 
point  with  linear-stability  theory. 

We  wish  to  emphasize  the  importance  of  knowing  the  proper 

value  of  As  noted  earlier,  uncertainty  about  the  value  of 

iiie  has  left  the  turbulence-model  transition-prediction  method 

with  an  effective  "adjustable  parameter"  in  previous  analyses. 

Removing  this  uncertainty  is  an  important  conceptual  advance. 

It  is  also  important  to  note  that  the  value  of  u)e  is  frequency 

dependent  so  that,  as  would  be  expected  on  physical  grounds, 

the  latter  stages  of  transition  are  predicted  to  be  not  entirely 

independent  of  spectral  effects.  While  we  have  yet  to  find  a 

7/4 

physical  explanation  for  the  propositionality  of  a)e  and  Fr  , 
we  continue  to  seek  such  an  explanation. 

In  conclusion,  results  to  date  show  great  promise  for  develop- 
ment of  a physically-sound  alternative  to  the  e9  procedure. 
Further  research  should  be  conducted  to  determine  how  well  the 
method  predicts  transition-point  location  for  well-documented 
flows  with  adverse  and/or  favorable  pressure  gradient.  Also, 
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further  testing  is  needed  to  more  definitively  establish  the 
method's  sensitivity  to  the  point  of  initiation  of  the 
turbulence-model  phase  of  the  computation.  When  such  analyses 
have  been  done  we  feel  confident  a reliable  analytical  tool 
for  predicting  boundary-layer  transition  will  result. 
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LIST  OF  SYMBOLS 


SYMBOL 


u,v,w 

u. 

Ue 

x,y,z 


8 , B* 


ir(y) 


DEFINITION 

Turbulent  mixing  energy  defined  in  Equation  (13) 
Dimensionless  streamf unction 
Dimensionless  frequency,  vfi/U^ 

Falkner-Skan  profile  parameter 
Static  pressure 
Plate-length  Reynolds  number 

Turbulent  Reynolds  number  defined  in  Equation  (8) 
Transition  Reynolds  number  based  on  momentum  thickness 
Reynolds  number  based  on  displacement  thickness 
Time 

Freestream  turbulence  intensity 
Velocity  components  in  x,y,z  directions 
Freestream  velocity 
Boundary-layer-edge  velocity 

Cartesian  coordinates  in  streamwise,  normal,  lateral 
directions 

Closure  coefficients 
Complex  wave  number 
Closure  coefficients 
Complex  wavenumber 
Falkner-Skan  profile  parameter 
Boundary-layer  thickness 
Kinematic  eddy  viscosity 
Blasius  similarity  variable 
Falkner-Skan  similarity  variable 
Closure  coefficient 

Average  value  of  X defined  in  Equation  (3D 

Minimum  value  of  X 

Kinematic  molecular  viscosity 

Complex  pressure  eigenfunction 

Density 

Closure  coefficients 
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LIST  OF  SYMBOLS 


SYMBOL 


DEFINITION 


<P(y) 

<t> 

U) 

Q 

l 


Complex  normal-velocity  eigenfunction 

Dimensionless  dissipation-rate  function  defined 
in  Equation  (33) 

Turbulent  dissipation  rate  defined  in  Equation  ( 14 ) 
Frequency 

Turbulent  length  scale,  e^/w 


Subscripts  and  Superscripts 

e Boundary-layer-edge  value 

i Imaginary  part 

r Real  part 

Other  Notation 

For  a given  variable  \Jr. 

Long-time-averaged  value  of  ^ defined  in  Equation  (20) 
Fluctuating  part  of  <p 


<\p> 
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